Insights into vertebral band pair deposition rate in the juvenile common thresher shark (Alopias vulpinus) in the northeastern Pacific Ocean.

Validation of band pair deposition rates in elasmobranch vertebrae is essential for accurate age estimation using band pair counting techniques. We present a validation study of the vertebral band pair deposition rate for juvenile common thresher sharks Alopias vulpinus in the northeastern Pacific Ocean (NEPO) using tag and recapture with oxytetracycline (OTC) injection. A total of 14 juvenile A. vulpinus marked with OTC from 1998 through 2013 were recaptured with times at liberty ranging from 1.08 to 3.81 years with an average of 2.14 years (±0.97 years standard deviation, SD). Shark size ranged from 80 to 128 cm fork length (LF) at the time of OTC injection and from 112 to 168 cm LF for those measured at recapture. The slopes of the relationships between band pairs post OTC and years at liberty for each reader ranged from 0.84 to 0.95, slightly lower than the 1.0 slope expected from annual band pair formation. These findings preliminarily support previous age and growth assumptions based on a one band pair per year deposition rate. However, high variation in band pair deposition rates between samples, coupled with regression slopes falling just under one band pair per year, indicates that further investigation is needed to refine band pair deposition rate estimates.

Natal origin of Pacific bluefin tuna from the California Current Large Marine Ecosystem.

Natal origin of subadult (age-1) Pacific bluefin tuna (PBT, Thunnus orientalis) from the California Current Large Marine Ecosystem (CCLME) was determined using natural tracers in ear stones (otoliths). Age-0 PBT collected from the two known spawning areas in the western Pacific Ocean (East China Sea, Sea of Japan) were used to establish baseline signatures from otolith cores over 4 years (2014-2017) based on a suite of trace elements (Li, Mg, Mn, Sr, Zn and Ba). Distinct chemical signatures existed in the otolith cores of age-0 PBT collected from the two spawning areas, with overall classification accuracy ranging 73-93% by year. Subadult PBT collected in the CCLME over the following 4 years (2015-2018) were then age-class matched to baselines using mixed-stock analysis. Natal origin of trans-Pacific migrants in the CCLME ranged 43-78% from the East China Sea and 22-57% from the Sea of Japan, highlighting the importance of both spawning areas for PBT in the CCLME. This study provides the first estimates on the natal origin of subadult PBT in this ecosystem using otolith chemistry and expands upon the application of these natural tracers for population connectivity studies for this species.

Species-Specific Characteristics Influence Contaminant Accumulation Trajectories and Signatures Across Ontogeny in Three Pelagic Shark Species.

Factors influencing organic contaminant accumulation in sharks, especially across ontogeny, are not well-known. Contaminant concentrations were measured in three species of sharks (Blue, Shortfin Mako, and Common Thresher) across a range of size classes (neonatal to adult) that vary in their ecological and physiological characteristics. Empirical data was compared to a theoretical framework that predicted the shape of lifetime accumulation curves. We found that a one-size-fits-all accumulation model was not appropriate as species-specific characteristics had a significant effect on contaminant accumulation trajectories. Maternal offloading likely has an important effect on determining neonatal shark contaminant starting points, and trophic ecology and physiology may interact to affect the shape of species' contaminant accumulation curves. Makos were found to have the highest accumulation potential and Blues the lowest, with Threshers being intermediate in accumulation potential. Changes in species' ecology and/or physiology were also reflected in contaminant signature changes over ontogeny. If contaminant concentrations are to be used as a proxy for risk, species-specific characteristics need to be taken into account when estimating contaminant exposure and its potential negative effects on shark health and human consumption safety.

Elements of time and place: manganese and barium in shark vertebrae reflect age and upwelling histories.

As upper-level predators, sharks are important for maintaining marine food web structure, but populations are threatened by fishery exploitation. Sustainable management of shark populations requires improved understanding of migration patterns and population demographics, which has traditionally been sought through physical and/or electronic tagging studies. The application of natural tags such as elemental variations in mineralized band pairs of elasmobranch vertebrae cartilage could also reveal endogenous and exogenous processes experienced by sharks throughout their life histories. Here, elemental profiles were characterized in vertebrae encompassing complete life histories (birth-to-death) of shortfin mako (Isurus oxyrinchus), common thresher (Alopias vulpinus) and blue shark (Prionace glauca) of known tag and recapture locations in the eastern North Pacific Ocean. All sharks were injected with oxytetracycline at initial capture, released and subsequently recaptured, with individual liberty times ranging from 215 days to 6 years. Vertebral band pairs forming over the liberty intervals were verified by counting the number of band pairs deposited since the oxytetracycline band. Regular oscillations in vertebrae manganese (Mn) content corresponded well with the number of validated band pairs, suggesting that Mn variation could be used to age sharks. Increases in vertebrae barium concentration were correlated with times when individuals occupied areas with high coastal upwelling indices, the timing and spatial intensity of which varied from year to year. Interspecific relationships were probably influenced by behavioural differences in horizontal and vertical habitat use, feeding habits and thermoregulatory physiology. These results indicate that vertebral sclerochronology has the potential to advance our knowledge of elasmobranch life history including age and growth estimation and environmental reconstruction.

The political biogeography of migratory marine predators.

During their migrations, marine predators experience varying levels of protection and face many threats as they travel through multiple countries' jurisdictions and across ocean basins. Some populations are declining rapidly. Contributing to such declines is a failure of some international agreements to ensure effective cooperation by the stakeholders responsible for managing species throughout their ranges, including in the high seas, a global commons. Here we use biologging data from marine predators to provide quantitative measures with great potential to inform local, national and international management efforts in the Pacific Ocean. We synthesized a large tracking data set to show how the movements and migratory phenology of 1,648 individuals representing 14 species-from leatherback turtles to white sharks-relate to the geopolitical boundaries of the Pacific Ocean throughout species' annual cycles. Cumulatively, these species visited 86% of Pacific Ocean countries and some spent three-quarters of their annual cycles in the high seas. With our results, we offer answers to questions posed when designing international strategies for managing migratory species.

An analytical approach to sparse telemetry data.

Horizontal behavior of highly migratory marine species is difficult to decipher because animals are wide-ranging, spend minimal time at the ocean surface, and utilize remote habitats. Satellite telemetry enables researchers to track individual movements, but population level inferences are rare due to data limitations that result from difficulty of capture and sporadic tag reporting. We introduce a Bayesian modeling framework to address population level questions with satellite telemetry data when data are sparse. We also outline an approach for identifying informative variables for use within the model. We tested our modeling approach using a large telemetry dataset for Shortfin Makos (Isurus oxyrinchus), which allowed us to assess the effects of various degrees of data paucity. First, a permuted Random Forest analysis is implemented to determine which variables are most informative. Next, a generalized additive mixed model is used to help define the relationship of each remaining variable with the response variable. Using jags and rjags for the analysis of Bayesian hierarchical models using Markov Chain Monte Carlo simulation, we then developed a movement model to generate parameter estimates for each of the variables of interest. By randomly reducing the tagging dataset by 25, 50, 75, and 90 percent and recalculating the parameter estimates, we demonstrate that the proposed Bayesian approach can be applied in data-limited situations. We also demonstrate how two commonly used linear mixed models with maximum likelihood estimation (MLE) can be similarly applied. Additionally, we simulate data from known parameter values to test each model's ability to recapture those values. Despite performing similarly, we advocate using the Bayesian over the MLE approach due to the ability for later studies to easily utilize results of past study to inform working models, and the ability to use prior knowledge via informed priors in systems where such information is available.

Fit to predict? Eco-informatics for predicting the catchability of a pelagic fish in near real time.

The ocean is a dynamic environment inhabited by a diverse array of highly migratory species, many of which are under direct exploitation in targeted fisheries. The timescales of variability in the marine realm coupled with the extreme mobility of ocean-wandering species such as tuna and billfish complicates fisheries management. Developing eco-informatics solutions that allow for near real-time prediction of the distributions of highly mobile marine species is an important step towards the maturation of dynamic ocean management and ecological forecasting. Using 25 yr (1990-2014) of NOAA fisheries' observer data from the California drift gillnet fishery, we model relative probability of occurrence (presence-absence) and catchability (total catch per gillnet set) of broadbill swordfish Xiphias gladius in the California Current System. Using freely available environmental data sets and open source software, we explore the physical drivers of regional swordfish distribution. Comparing models built upon remotely sensed data sets with those built upon a data-assimilative configuration of the Regional Ocean Modelling System (ROMS), we explore trade-offs in model construction, and address how physical data can affect predictive performance and operational capacity. Swordfish catchability was found to be highest in deeper waters (>1,500 m) with surface temperatures in the 14-20°C range, isothermal layer depth (ILD) of 20-40 m, positive sea surface height (SSH) anomalies, and during the new moon (<20% lunar illumination). We observed a greater influence of mesoscale variability (SSH, wind speed, isothermal layer depth, eddy kinetic energy) in driving swordfish catchability (total catch) than was evident in predicting the relative probability of presence (presence-absence), confirming the utility of generating spatiotemporally dynamic predictions. Data-assimilative ROMS circumvent the limitations of satellite remote sensing in providing physical data fields for species distribution models (e.g., cloud cover, variable resolution, subsurface data), and facilitate broad-scale prediction of dynamic species distributions in near real time.

Prediction of fishing effort distributions using boosted regression trees.

Concerns about bycatch of protected species have become a dominant factor shaping fisheries management. However, efforts to mitigate bycatch are often hindered by a lack of data on the distributions of fishing effort and protected species. One approach to overcoming this problem has been to overlay the distribution of past fishing effort with known locations of protected species, often obtained through satellite telemetry and occurrence data, to identify potential bycatch hotspots. This approach, however, generates static bycatch risk maps, calling into question their ability to forecast into the future, particularly when dealing with spatiotemporally dynamic fisheries and highly migratory bycatch species. In this study, we use boosted regression trees to model the spatiotemporal distribution of fishing effort for two distinct fisheries in the North Pacific Ocean, the albacore (Thunnus alalunga) troll fishery and the California drift gillnet fishery that targets swordfish (Xiphias gladius). Our results suggest that it is possible to accurately predict fishing effort using < 10 readily available predictor variables (cross-validated correlations between model predictions and observed data -0.6). Although the two fisheries are quite different in their gears and fishing areas, their respective models had high predictive ability, even when input data sets were restricted to a fraction of the full time series. The implications for conservation and management are encouraging: Across a range of target species, fishing methods, and spatial scales, even a relatively short time series of fisheries data may suffice to accurately predict the location of fishing effort into the future. In combination with species distribution modeling of bycatch species, this approach holds promise as a mitigation tool when observer data are limited. Even in data-rich regions, modeling fishing effort and bycatch may provide more accurate estimates of bycatch risk than partial observer coverage for fisheries and bycatch species that are heavily influenced by dynamic oceanographic conditions.

No effect of trans sodium crocetinate on maximal O(2) conductance or V(O(2),max) in moderate hypoxia.

The lumped parameter describing skeletal muscle diffusional conductance for O(2), DM(O(2)), reflects all of the resistances for O(2) in moving from red cell to muscle fiber mitochondria. The purpose of our study was to determine if the carotenoid compound, trans sodium crocetinate (TSC), which has been reported to increase the diffusivity of O(2) in plasma, improves DM(O(2)) and thus, V(O(2),max) in maximally contracting in situ skeletal muscle. V(O(2),max) was measured in the isolated perfused canine gastrocnemius (n=5) during 3 min of isometric tetanic contractions at 1 Hz, while the animal was breathing 12% O(2) (PA(O(2))=32+/-2 Torr, mean+/-S.E.) under two experimental conditions. The first was a control contraction period and the second (following 60 min recovery) was performed within 5 min after infusion of a 0.1 mg x ml(-1) solution of TSC (total dose 100 microg kg(-1)). There were no significant differences in convective O(2) delivery (11.9+/-2.3 vs. 12.1+/-2.2 ml x min(-1) x 100 g(-1)), V(O(2),max) (9.5+/-1.5 vs. 9.6+/-1.5 ml x min(-1) x 100 g(-1)) or calculated DM(O(2)) (0.37+/-0.03 vs. 0.37+/-0.04 ml x min(-1) x 100 g(-1) x Torr(-1)) between contraction periods. As such, our results show that TSC does not improve performance in maximally contracting canine gastrocnemius muscle in situ under moderately hypoxic conditions, suggesting either that TSC in this situation does not increase plasma O(2) diffusivity or that this step in O(2) diffusion from red cell to myocyte does not constrain DM(O(2)).